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Abstract 

Gribov theory is applied to investigate the shadowing effects in the structure functions 
of nuclei. In this approach these effects are related to the process of diffractive dissociation 
of a virtual photon. A model for this diffractive process, which describes well the HERA 
data, is used to calculate the shadowing in nuclear structure functions. A reasonable 
description of the x, and A-dependence of nuclear shadowing is achieved. 
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1. Introduction 

Deep inelastic scattering (DIS) on nuclei gives important information on distributions 
of quarks and gluons in nuclei. The region of small Bjorken x is especially interesting 
because partonic clouds of different nucleons overlap as a; — > and shadowing effects 
become important. There are experimental results in this region, which show that there 
are strong deviations from an behavior in the structure functions [1] . Several theoretical 
models have been proposed to understand these data [1]. The most general approach is 
based on the Gribov theory [2]. It relates partonic and hadronic descriptions of small 
X phenomena in interactions of real or virtual photons with nuclei. In this approach the 
shadowing effects can be expressed in terms of the cross-sections for diffraction dissociation 
of a photon on a nucleon (Fig. 1). This process has been studied recently in DIS at 
HERA [3]. The detailed and (M is the invariant mass of the diffractively 

produced system) dependencies observed in these experiments have been well described 
in the theoretical model of ref. [4] which is based on Regge factorizations and uses as 
an input available information on diffr active production in hadronic interactions. Here 
we will apply the same model to calculate the structure functions of nuclei in the small 
x-region. The use of the model, which describes well the diffraction dissociation of virtual 
photons on a nucleon target, leads to a strong reduction of the theoretical uncertainty in 
calculations of the structure functions of nuclei in comparison with previous calculations 
[1, 5-8]. It also allows to discuss the shadowing effects in gluon distributions. 

2. The model 

In the Gribov approach the forward scattering amplitude of a photon with virtuality 
on a nuclear target can be written as the sum of the diagrams shown in Fig. 2. Since 
we are interested in the low x region we will describe the various interactions by 



Pomeron exchange. The diagram of Fig. 2a corresponds to the sum of interactions with 
individual nucleons and is propotional to A^. The second diagram (2b) contains a double 
scattering with two target nucleons. It gives a negative contribution to the total cross- 
section, proportional to A"^/^ (for large A). It describes the first shadowing correction 
for sea quarks. According to reggeon diagram technique [9] and Abramovsky, Gribov, 
Kancheli (AGK) cutting rules [10], the contribution of the diagram of Fig. 2b to the total 
^*A cross-section is related to the diffractive production of hadrons by a virtual photon 
as follows : 
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where Ta(&) is the nuclear profile function, pA is the nuclear density (7^(6) = dZpA{b, Z), 
JdHTAib) = 1) and 

FA{tmin) = j d^hJQ{^\/ -tminb)TA(h) , tmin = -m% (^ M^+ ) 

Note that FA{tmin) is equal to unity as X — > and decreases fast as x increases to x^^ 
^ ^r, , due to a lack of coherence for x > Xcr- 

Eq. (1) is written in the approximation R\ ^ Rj^, where Rn is the radius of the 
7*p interaction. It will be used in this form only for ^ > 20 (see below). We have 
also neglected the real part of the Pomeron amplitude which is small for our value of 
the Pomeron intercept (see Eq. (5)). However, for higher values of this intercept the 
contribution of the real part can be substantial [11]. 

For a deuteron, the double rescattering contribution has the following form 
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where FD{t) — exp(at), with o = 40 GeV~^. M^^^ in eqs. (1), (2) corresponds to the 
minimal mass of the diffr actively produced hadronic system and M^^^, is chosen according 
to the condition : xp = x- ^'j",^' ~ ^^J^' < 0.1. 



Equation (2) has been used to calculate inelastic contributions to Glauber corrections 
in hadron-deuteron interactions [12, 13] and was generalized to heavier nuclei in the form 
(1) in ref. [14]. 

Thus the second order rescattering term can be calculated if the differential cross- 
section for diffractive production by a virtual photon is known. 

Higher order rescatterings are model dependent, but calculation shows that, for the 
values of A and x {x > 10~^) where experimental data exist, their contribution is rather 
small. We use the following unitary expression for the total cross-section 



where 



dM'^dt 

t=o 

This expression is valid in the generalized Schwimmer model [15, 16] and is obtained from 
a summation of fan diagrams with triple Pomeron interaction. However, its physical basis 
and applicability is much broader. For example it follows from the rescattering of a qq 
system with transverse sizes distributed according to a gaussian [17]. We have checked 
that the results obtained from a summation of higher order rescatterings of an eikonal 
type are very similar to the ones obtained with (3) - the differences being of the order of 
one percent. 

Thus we have for the ratio Ra = F2A/F2N of nucleus and nucleon structure functions, 
in the region of small x 

F2A _ [ 2. ATAjb) 

F2N J l+{A-l)f{x,Q^)TA{b) • 

The deviation of this ratio from A^ = A J d'^hTA{h) is due to the second term in the de- 
nominator of the integrand in eq. (4). Thus, knowing the differential cross-section for 
diffraction dissociation on a nucleon and the structure function of a nucleon (cr^*jv), one 
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can predict the A (and x, Q^) dependence of structure functions of nuclei. Eq. (4) can only 
be used in the region of x substantially smaller than 10~^ where the sea quarks component 
dominates. For x close to lO""*^ shadowing of valence quarks (which in general is not de- 
scribed by eq. (4)) becomes important [18, 19]. The effects leading to antishadowing (such 
as real parts in the rescattering diagram due to secondary exchanges) are also important 
in the region of a; ~ 0.1. 

In refs. [4] we described the diffractive contribution to DIS in terms of Pomeron 
exchange 

F^ix, Q^ xp, t) = ^P^x'f'^-^'^FpiP, Q\ t) (5) 

where gpp{t) is the Pomeron-proton coupling {gpp{t) = 5f^(0) exp(Ct) with {gpp{0))^ = 
23 mb and C = 2.2 GeV-^), ap{t) = ap(0) + a'p(0)t is the Pomeron trajectory {ap{0) 
= 1.13, a'p{0) = 0.25 GeV-2) and Fp{(3,Q'^,t) is the Pomeron structure function. The 
variable /? = m^^^q"^ ~ ^ plays the same role for the Pomeron as the Bjorken variable x 
for the proton. At large Q^, Fp can be expressed in terms of the quark distributions in 
the Pomeron 

Fp(/3, t) = J2 4P [qfiP, Q\t) + t)] ■ (6) 

i 

In refs. [4] we determined Fp{P, Q^, t) using Regge- factorization for small values of P and 
a plausible assumption on the /5 — > 1 behavior. This function was then used as an initial 
condition for QCD evolution of partons in the Pomeron. The results of the QCD-evolution 
crucially depend on the form of the gluon distribution in the Pomeron. Experimental 
results for F2' can be understood only if the distribution of gluons in the Pomeron is 
rather hard and the gluons carry the main part of the Pomeron momentum [4, 20-22] . The 
explicit forms of all these functions are given in Appendix 1. 

The validity of Pomeron factorization (5) for F^ as well as that of the QCD evolution 
for partons in the Pomeron has been questioned in recent papers. These papers deal with 
diffractive charm production [23, 24] and with the contribution of longitudinal photons to 
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difFractive production [25, 26]. However, in all these papers high-twist effects (in or Mq, 



where Mq is the mass of the heavy quark), which give small contributions to difFractive 
cross-sections, were considered. Arguments in favour of usual QCD evolution for the main 
twist contribution to F^" have been given in ref. [27]. In any case the CKMT model [4] 
gives a reasonable description of difFractive production in DIS. Thus it effectively includes 
high twist effects and can be used to compute the function f{x, Q^), which determines the 
shadowing of nuclear structure function via (4). This function can be written in terms of 
the ratio Fp/F2n '■ 



where the integration limits are x/xqp with xqp = 0.1 and Q'^/{M^^^ + Q"^). In the 
following we take M^^^ = 0.4 GeV^, in order to include the p- meson peak in the integration 
region. 

The parametrization of the Pomeron [4] and nucleon [28] structure functions are given 
in Appendix 1. Note that the Q^-dependence of nuclear shadowing is obtained by evolving 
separately the nucleon and Pomeron structure functions and taking their ratio in eq. (7). 
Actually, one should compute first F2A at = Qq and evolve it using the nuclear partonic 
distributions. However, this would require the knowledge of these distributions for all 
values of x. At small x, where sea quarks are dominant, these two procedures are equivalent 
for the Born term and the first rescattering correction in eq. (4) . As discussed above, higher 
rescattering corrections are small. 

In the numerical calculations we use a standard Woods-Saxon profile Ta(&) for A > 20. 
For light nuclei {A < 20) we use a gaussian profile 





TA{b) 



3 



(8) 



27tR\ 



with an r.m.s. radius parametrized as [29] 
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Ra = 0.82 A^/^ + 0.58 fm 



(9) 



For deuteron eq. (9) is not valid. In this case we use eq. (2). The simple exponential 
form of -Fd(^) gives results which differ by less than 20 % from the more sophisticated 
parametrization used in refs. [6] [30] . 

In eq. (1) we have neglected the t-dependence of the 7*p diffractive cross-section. As 
explained above, this approximation is only used for large nuclei where nucleon sizes can 
be neglected as compared to nuclear ones. For light nuclei {A < 20), we take into account 
this t (or 6)-dependence by making the following replacement 



This nucleon radius approximately describes the t-dependence of the 7*p diffractive cross- 
section in the kinematical region we are interested in. 

3. Numerical results 

The results of our calculations are shown in Figs. 3-8. Theoretical predictions for the 
deuteron structure function /2F.^ are shown in Fig. 5. Our results are close to those 
of refs. [30, 31] but smaller by a factor of about 3 from the results of ref. [6]. Comparison 
of our predictions for the ratio \F2^ /F^ with experimental data of NMC [32] is shown in 
Fig. 3 and for ratios of different nuclei in Fig. 4. New data for the ratio F^'^ jF^ [33] are 
also shown in Fig. 4. It is important to note that experimental points in Figs. 3, 4 for 
different x correspond to different values of < > [32] [33]. This correlation has been 
taken into account in our calculations. The agreement between theoretical predictions and 
experimental data is good. Note that there are no free parameters in our calculations. 



= 0.8 fm 



(10) 
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Our predictions for the ratio ^ in the region of very small x are shown in Figs. 5 
for fixed values of Q^. They could be confronted to experiment if nuclei at HERA would 
be available. Note that our results are more reliable for small values of a; (a; < 10~^) where 
the effects of both valence quark shadowing and antishadowing are negligeable. The curves 
for shadowing effects in the gluon distribution of nuclei ^^w^^^r^ are shown in Figs. 6. 
They look similar to the shadowing for the quark case. However the absolute magnitude 
of the shadowing is smaller in the gluon case contrary to expectations of some theoretical 
models [1] but in agreement with [34] . (Note that these results are sensitive to the gluonic 
distribution in the Pomeron, which is poorly known at present). These predictions can be 
tested in experimental studies of J/ip and T-production on nuclear targets at RHIC and 
LHC. 

Finally we want to discuss in more detail the Q^-dependence of the shadowing. Recent 
NMC data [33] for the ratio of F^"/F^ are shown in Figs. 7 as functions of for fixed 
values of x in the small x region. The theoretical curves have a weak dependence on 
and are in reasonable agreement with experiment, although the dependence seems 
stronger in the data especially in the region of small Q^. At larger values of x the data are 
practically Q^-independent. These properties should be checked in future experiments. 

4. Conclusions 

In conclusion, a model based on the Gribov-Glauber theory of nuclear shadowing 
and the properties of diffraction in DIS observed at HERA, leads to a fair description of 
experimental data on structure functions of nuclei in the small x region. Predictions of 
shadowing effects for quark and gluon distributions are given. They can be tested in future 
experiments at HERA and in hadronic colliders. 



9 

Acknowledgement s 

One of us (A.C.) would like to thank G. Do Dang for discussions. A. K. and D. P. wish 
to thank the LPTHE for hospitality during a period when this work was initiated. This 
work has been partially supported by grant 93-0079 of INTAS. A. K. also acknowledges 
support from the grant N° 96-02-19184 of RFFI. 



Appendix 1 



10 



For the Pomeron and nucleon structure functions we use the parametrization of the 
CKMT model [4, 28] 

F2n{x, = yl(g2)3,-A(Q^)(i _ 3,)n(Q^)+4 ^ ^^^^ Q^)x^-<-^^^\l - a;)"^^') , {A.l) 

Fp{P, Q2) = F2n{P, g'; A ^ eA, B{x) ^ fB', n^n-2) {A.2) 

where 

with (all dimensional quantities are in Gev^) 

^ = 0.1502 , S' = 1.2035 , = 0.4150 , Aq = 0.07684 

a = 0.2631 , 6 = 0.6452 c = 3.5489 , d = 1.1170 , e = / = 0.07 . 
Finally, we have [28] 

S(a;) = 0.754 + 0.4495(1 -x) . (A.3) 

The two terms in Eq. (A. 3) appear because we have used for the valence quark distributions 
in the proton d{x) = u{x){l — x). Such a difference between u and d quark distributions 
is not present in the Pomeron case, and we have dropped the 1 — x factor in (A. 3). In the 
CKMT model the Pomeron structure function is determined from the nucleon one using 
triple Regge couplings determined from soft diffraction data and Regge factorization. 

Comparison of eqs. (A.2) and (6) allows to determine the valence and sea quark 
distributions in the Pomeron. Likewise one can determine the corresponding ones in the 
nucleon. These distributions are used as initial conditions at = Qq in the DGLAP 
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evolution equation, as described in [28]. The gluon distributions for < Qq in the 
nucleon and the Pomeron are [4, 28] 

gAx,Q^) = e AgiQ^) x-^^'^'^l - x)-'-' {A.b) 
where Ag[Q'^) has the same form as A[Q'^) in eqs. (A.l) and (A. 2) : 

The normahzation of is obtained from the energy-momentum sum-rule. For the 
Pomeron this sum-rule is not valid and the normalization of gp is obtained from that of 
g^ using Regge factorization. The constant e = 0.07 is the same as in eq. (A. 2). Actually 
there is a large ambiguity in the shape of gp. We have only determined the x-behaviour at 
small X as well as the absolute normalization. The form (A. 5) is just a simple extrapolation 
to the region of a; ^ 1. 

Using the gluon distributions in eqs. (4) and (7) one can obtain the corresponding 
distributions in nuclei. 

In order to have exactly the diffractive cross-section computed in ref. [4] as well as 
the same F2N of ref. [28] we use the values of Qq in those references. These are Qq = 
5 GeV^ for the Pomeron and Qq = 2 GeV^ for the nucleon. The corresponding gluon 
normalizations obtained from the energy- momentum sum rule are Ag — 1.84 at Qq = 
2 GeV2 and Ag = 1.71 at = 5 GeV^. 



12 

References 

[1] M. Arneodo, Phys. Reports 240 (1994) 301 (and references therein). 

[2] V. N. Gribov, ZhETF 56 (1969) 892, ibid 57 (1969) 1306 [Sov. Phys. JETP 29 
(1969) 483, 30 (1970) 709]. 

[3] T. Ahmed et al (HI coUaboration), Phys. Lett. B348 (1995) 681. 

M. Derrick et al (Zeus collaboration), Z. Phys. C68 (1995) 569 ; Z. Phys. C70 (1996) 
391. 

[4] A. CapeUa, A. Kaidalov, C. Merino and J. Tran Thanh Van, Phys. Lett. B343 (1995) 
403. 

A. Capella, A. Kaidalov, C. Merino, D. Pertermann and J. Tran Thanh Van, Phys. 
Rev. D53 (1996) 2309. 

[5] K. Boreskov, A. CapeUa, A. Kaidalov and J. Tran Thanh Van, Phys. Rev. D47 
(1993) 219. 

[6] V. Barone et al, Z. Phys. C58 (1993) 541. 

[7] W. Melnitchouk, A. W. Thomas, Phys. lett. B317 (1993) 437 ; Phys. Rev. C52 
(1995) 3373. 

[8] G. Piller, W. Ratzka, W. Weise, Z. Phys. A352 (1995) 427. 

[9] V. N. Gribov, ZhETF 57 (1967) 654 [Sov. Phys. JETP 26 (1968) 14]. 

[10] V. A. Abramovsky, V. N. Gribov and O. V. Kancheh, Yad. Fiz. 18 (1973) 595 [Sov. 
J. Nucl. Phys. 18 (1974) 308]. 

[11] A. Bialas, W. Czyz and W. Florkowski, TPJU-25/96 and references therein. 

[12] A. B. Kaidalov, L. A. Kondratyuk, JETP Letters 15 (1972) 170 ; Nucl. Phys. B56 
(1973) 90. 



13 

[13] V. V. Anisovich, L. G. Dakhno and P. E. Volkovitsky, Yad. Fiz. 15 (1972) 168 [Sov. 
J. Nucl. Phys. 15 (1972) 97]. 

[14] V. A. Karmanov and L. A. Kondratyuk, JETP Letters 18 (1973) 266. 

[15] A. Schwimmer, Nucl. Phys. B94 (1975) 445. 

[16] K. G. Boreskov et al., Yad. Fiz. 53 (1991) 569 [Sov. J. Nucl. Phys. 53 (1991) 356]. 
[17] B. Kopeliovich et al., JINR E2-86-125. 

[18] L. L. Frankfurt, M. I. Strikman and S. Liuti, Phys. Rev. Lett. 65 (1990) 1725. 

[19] A. B. Kaidalov, C. Rasinariu and U. Sukhatme, UICHEP-TH/96-9. 

[20] T. Gehrmann and W. J. Stirling, Z. Phys. C70 (1996) 89. 

[21] K. Golec-Biernat and J. Kwiecinski, Phys. Lett. B353 (1995) 329. 

[22] J. Dainton (HI collaboration). Proceedings Workshop on Deep Inelastic Scattering 
and QCD, Paris, France, 24-28 April 1995 (ed. J. P. Laporte and Y Sirois). 

[23] M. Genovese, N. N. Nikolaev, B. G. Zakharov, Phys. Lett. B378 (1996) 347. 

[24] E. M. Levin, A. D. Martin, M. G. Ryskin, T. Teubner, DTP 96-50. 

[25] J. Bartels, H. Letter, M. Wusthoff, Phys. Lett. B379 (1996) 248 ; J. Bartels, C. 
Ewerz, H. Lotter, M. Wusthoff, Phys. Lett. B386 (1996) 389. 

[26] M. Genovese, N. N. Nikolaev, B. G. Zakharov, Phys. Lett. B380 (1996) 213. 

[27] A. Berera and D. E. Soper, Phys. Rev. D53 (1996) 6162. 

[28] A. CapeUa, A. Kaidalov, C. Merino, J. Tran Thanh Van, Phys. Lett B337 (1994) 
358. 

[29] M. A. Preston and R. K. Bhoduri, Structure of the Nucleus, Addison- Wesley, New 
York 1975. 

[30] W. Melnitchoak, A. W. Thomas, Phys. Rev. D47 (1993) 3783. 



14 

[31] B. Badelek, J. Kwiecinski, Phys. Rev. D50 (1994) 4. 
[32] P. Amandruz et al (NMC collaboration), Nucl. Phys. B441 (1995) 3. 
[33] N. Arneodo et al (NMC collaboration), Nucl. Phys. B481 (1996) 23. 
[34] K. J. Eskola, Nucl. Phys. B400 (1994) 240. 



Figure Captions 
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Fig. 1 : DiflFractive dissociation of a virtual photon. The shaded area represents the 
exchange of a Pomeron. 

Fig. 2 : The first two terms (single and double scattering) of the multiple scattering 
series for the total 7*A'" cross-section in the Gribov-Glauber approach. 

Fig. 3 : The ratios (^IA)F^lF^ computed from eq. (3) for different values of x. The 
experimental points are from ref. [17]. The values of are different for 
different^- values [17]. 

Fig. 4 : The ratios {^Axj A-2)F^'^ j F^^ computed from eq. (3) for different values of x. 

The experimental points are from refs. [17] and [18]. The values of are 
different for different x values [17, 18]. 

Fig. 5 : The ratios [Yj A)F^ jF^ computed from eq. (3) for different values of x in the 
small X region, at fixed values of Q^. 

Fig. 6 : The ratios [Xj A'i)g^ j of gluon distribution functions computed for different 
values of x in the low x region, at fixed values of . 

Fig. 7 : The ratio (12/119)F2^"/F2^ computed from eq. (3) for different values of Q^, 
at two fixed values of x. The data points are from ref. [18]. 
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